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The details of the calculation of the two-boson exchange effects in the parity- 
violating elastic ep scattering within a simple hadronic model, including both the 
nucleon and A(1232)-resonance intermediate states, are presented. We examine the 

m 

' sensitivity of our results with respect to choice of form factors. We emphasize the 

' importance to use correct relations relating — )• A and A ^ N transition vertex 

functions. The A^A Coulomb quadrupole transition is found to play important role 
at higher > 3.0 GeV^. We also elucidate the relation between our results and the 
' well-known result on the effect given by Marciano and Sirlin (MS). The effect 



of the nucleon contribution 5n to parity-asymmetry Apy, is found to be in general, 
larger than the corresponding A contribution 5/\ except at extreme forward angles. 
The corrections to the extracted values of the strange form factors -|- /3G^^ from 
the HAPPEX, A4, and GO data are also presented. The total TBE corrections to 
the extracted values of -|- /3G^ in recent experiments of HAPPEX GO, and A4 
are, depending on kinematics, found to be small except in a few cases where they 
range from -20.6% to 48.3%. 
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I. INTRODUCTION 



One of the most intriguing questions in hadron structure is the possible existence of 
strangeness content in the proton since practically all constituent quark models employ only 
u and d quarks for light baryons. It was prompted by the EMC experiments l| which 
indicate that the amount of spin carried by the strange quark pairs ss is comparable to that 
carried by the u and d quarks and polarized opposite to the nucleon spin. Similar conclusion 
was also drawn from elastic up scattering [2;j| and theoretical analysis of vrA^ sigma term j^. 
A few other experiments have since been proposed including the excess of production 
in pp annihilation jsj, A polarization in deep-inelastic neutrino scattering jo, 7|, and double 
polarizations in photo- and electroproduction of meson 8| scheduled at SPringS for 2010 
9|, and the parity- violating electron-proton scattering. 

Parity-violating ep scatterings was first suggested as a unique probe to extract proton 



strange form factors by Kaplan and Manohar 



10| from measuring the parity-violating asym- 



metry Apv = {ctr — ctl)/ {o-r + (Xl) with polarized electrons, where ctr^l) is the cross section 
with a right-handed (left-handed) electron. The asymmetry arises from the interference of 
weak and electromagnetic amplitudes. Weak neutral current elastic scattering is mediated 
by the Z-boson exchange and measures form factors which are sensitive to a different linear 
combination of the three light quark distributions. When combined with proton and neu- 
tron electromagnetic form factors and with the use of charge symmetry, the strange electric 
and magnetic form factors, G% and G%f, can then be determined jlO|. Since this is a rather 
clean technique to access the charge and magnetization distributions of the strange quark 



within nucleons, four experimental programs SAMPLE 



GO 



11 , HAPPEX 



12|, A4 (l3|, and 



14| have been designed to measure this important quantity, which is small and ranges 



from 1 to 100 ppm. These experiments have been able to reach a precision of 6Apy ~ 0.1 



ppm. Several global analyses have been performed |15l-ll7l| and found that the electric and 
magnetic strange form factors are quite small with considerable error bars. Accordingly, 
greater effort to reduce theoretical uncertainty is needed in order to arrive at a more reliable 
interpretation of experiments. 

Leading order radiative corrections to Apy, including the box diagrams Fig. 1(d) and 
other diagrams, have been extensively studied 18N21| and widely used in the global analyses 



m 



15- 



ITj . Among those corrections, the interference between 7Z exchange {'jZE) of Fig. 
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1(d) with Fig. 1(a), was evaluated within the zero momentum transfer approximation, i.e.. 



= 0. The first calculation beyond the = approximation was done in 22| where 
the contribution of the interference of the two-photon exchange {2'~fE) process of Fig. 1(c) 
with diagram of Figs. 1(a) and 1(b) to Apy, was evaluated in a partonic approach using 
GPDs. It was prompted by the fact that such a parton model calculation of the 2'jE effect 



23| was arguably able to quantitatively resolve the discrepancy between the measurements 



of the proton electric to magnetic form factor ratio R = ^pGe/Gm, where /ip = 2.79, 
from Rosenbluth technique and polarization transfer technique at high momentum-transfer- 
squared Im]. It was found 22 1 that the correction to Apy is both Q"^ and e 
dependent, and can reach several percent in certain kinematics, becoming comparable in 
size with existing experimental measurements of strange-quark effects in the proton neutral 



22 



23| are reliable on 
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y for large 
141 have been 



weak current. However, the partonic calculations of 
comparable to a typical hadronic scale, while all current experiments 
performed at lower values. 

The two-boson exchange (TBE) corrections to Apy, namely, the contributions of the 
interference of the two-photon exchange {2'~fE) process of Fig. 1(c) with diagram of Figs. 
1(a) and 1(b) to Apy, and that between the 'jZ exchange of Fig. 1(d) with Fig. 1(a), 
were investigated in a hadronic model first with only intermediate states restricted to elastic 



nucleon states in 



25 



26| . This hadronic model was developed in |27(] to evaluate the 2^E 



contribution to the ratio R. The advantage of such a hadronic approach [27[ is that it is 
applicable to low region and the results obtained are in agreement with the partonic 



26| that both the the 2'yE and jZE corrections to Apy 



calculation of [23|]. It is found 25 . 
depends strongly on and e, and can reach a few percent and are comparable in size with 



the current experimental measurements of strange quark effects in the proton weak neutral 
current and their combined effects on the extracted values of G|; + f3Gpj can be as large as 

26] that the results show some sensitivity 



—40% in certain kinematics. It was further found 
on whether a monopole or dipole form is assumed for the nucleon form factors. 

Recently, the hadronic calculations on the TBE effects 25|, l26|] were extended to include 



A(1232) resonance in the intermediate states 28|, |29|] since A(1232) is known to play a 



dominant role in low-energy hadron physics 30]. Both calculations show that the interplay 
between the nucleon and A contributions depend strongly on the kinematics. However, 
there are discrepancies in the size of the total TBE corrections due to the use of different 
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FIG. 1: (a) One-photon exchange, (b) Z-boson exchange, (c) Two-photon exchange, and (d) 7Z- 
exchange diagrams for elastic ep scattering. Corresponding cross-box diagrams are imphed. 

vertex relation relating the vertices of — )• A and 7A — )■ A^, the strength of the Coulomb 
quardrupole excitation of the A, and the A form factors. 



In this paper, we give the details of our hadronic model calculations [25|, 128[ of the 2^E 
and '-yZE corrections to Apy and present a more extensive results of our calculation. In 
)articular we analyze in details the difference between our calculations and those of Ref. 



26 



29| . In addition, we demonstrate explicitly that our results do recover the results of [2 
in the limit of = 0. 

This article is organized as follows. The formalism for parity-violating electron-proton 
is given in Section II. The details of our calculation of the '^ZE and 2'-)E box diagrams in 
a simple hadronic model are presented in Section III. The numerical results of the above 
calculations and the impacts of our results on the extraction of the strange form factors and 
the weak charge of the proton are discussed in Section IV. In section V, we summarize our 
work. 



II. PARITY- VIOLATING ELECTRON-PROTON ELASTIC SCATTERING 

In this section we first briefly present the formulation of the parity-violating electron- 
proton elastic scattering within one-boson exchange (OBE) approximation and the corre- 
sponding procedure to extract the proton strange form factors. We then go beyond the OBE 
framework and discuss radiative corrections. 
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A. Parity-violating ep scattering within one-boson-exchange approximation 



The OBE diagrams of the elastic electron-proton scattering, e{pi)+p{p2) ^(pa) +^(^4), 
include one-photon exchange (I7-E') and one-Z-boson exchange (IZE) diagrams, as shown in 
Fig. 1(a) and 1(b), respectively. At hadron level, the couplings of the photon and Z-boson 
with the proton are given as 



{p'\J'\p) = W) 

{pvr\p) = w) 



2M, 



N 



2Mn 



q 



u{p), 



(1) 



where is the proton mass and q = p' — p. and are the proton electromag- 

netic/neutral weak current and axial form factors, respectively. The Sachs form factors are 
defined as 

^ p,/z _ ^p,/z^ ^,/z ^ p,/z ^ p,/z^ (2) 
with = —q^. The OBE diagrams of Figs. 1(a) and 1(b) are given in 



where r 



terms of the matrix elements of the electromagnetic and neutral weak currents 



-iu{p^){-ie'^^)u{pi 



le 



-u{Pi)Tl{pi,P2)u{p2) 



-iu{p3 



4 cos 9 



w 



[(-l + 4sin2 ew)+l5Hpi) 



—I 



q^-Ml 



IE 



u{p4)'r^{p4,P2)u{p2) 



(3) 



where ri{p',p) = ze{N{p')\r;^\N{p)) , T^^{p',p) = {zg / 4 cos 9w){N{p')\J^\N{p)). g = 
e/ sin6]Y is the weak coupling constant with 6^/ the Weinberg weak mixing angle, and 
Mz the Z-boson mass. The parity asymmetry in OBE approximation arises from the inter- 
ference of M'-"^ and M*-^-*. Straightforward calculation, leads to the following expression of 
parity-asymmetry in OBE approximation in terms of the form factors defined in Eq. ([1]) 

GfQ^ Ae + Am + Aa 



(4) 



Aa = -{1-A sin^ Ow) ^T{l + T){l-e^)GlGl, , 

where e = [1 + 2(1 + r) tan^ 6'iafe/2]^^ and Oiab the scattering angle of the electron in the 
laboratory frame. Gp = V2g'^/8M^ = 7raem/(y2M| cos^ 9w sin^ Ow) = 1.166 x IQ-^GeV-^ 
is the Fermi constant and ctem = e^/47r the fine structure constant. 
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To extract the strange form factors from Eq. (jlj) , one needs to make flavor decompositions 
of the form factors and G^\j. In the standard model, the electromagnetic current 

and the neutral weak current are given as 



where 



3l|, 



jem 



f=u,i 



f=U,d,! 



(5) 



9'v 



-l + 4sin2^p^, (7i = +l, 



9v 



9'v 



9v 



1 --sin^Ow, 9a = -1, 
4 

-l + -sm''ew, 9A = +h 
4 

-1 + -sin^^w', 9l = -1- 



(6) 



From Eqs. ( ]T]l5]l6l) . one obtains 



9 1 
'^^E,M ~ '^^E,M 



'E,Mi 



Iqs/p 

3 

r^^E,M^ 



G 



Z,p 



yu/p 



EM ~ ^ 3 ^w)Ge'm 



4 
3 



-1 + ^ sin^ ew)Gf^ + + l sin^ Ow)G'^% 



(7) 



where G^.m defined as follows 



{p{p')\qnnqf\p{p)) = Mp') 



2M 



Up{p), 



(8) 



and 



(9) 



If charge symmetry is assumed, i.e., the distribution of the u quarks in the proton is the 
same as that of the d quarks in the neutron, then one has G^'^^j = G^'''^,^, G^^^^ = G^^.j, and 



G 



s/p 
E,M 



G 



s/n 
E,M 



G|; ^.]- such that we can express the neutron electromagnetic form factors 



as 



^7,n _ ^d/p _ ^ ^u/p _ 

^E,M — g'-^E.M r^^E,M r^^E,M- 



(10) 



Combining the first equation of Eq. ([7]) and Eq. (fTOj) leads to the following two relations 



G 



u/p 
EM 



-^^EM "T ^ EM 



^EM' ^EM 



^ EM 



'^^'em + ^ 



EM- 



(11) 
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Putting the above two relations of Eq. (fTTj) back to the last relation in Eq. ([7]), the neutral 
weak form factors can be expressed in terms of the electromagnetic form factors of the proton 
and neutron, and the strange form factors 



4sin2 ey,)Gl\,-G^^'^ 



EM 



G 



E,M' 



(12) 



With Eq. ( IT2ll . the parity asymmetry Apy of Eq. (14]) can be rewritten as, 



Apv 
A^ = 



Ai + A2 



eiGYY + riGiry 



' M 



a(l - 4sin^ 9w] 



r{Glfr 



(13) 



where a = G pQ'^ / 4:7Taem\^ and e' = a/t(1 + t)(1 — e^). The electromagnetic form factors 
G'em ^E^M "^^^ t)e extracted from the elastic electron scattering from proton and 
deuteron (for the neutron), and the axial form factor can be extracted from the pion 



photoproduction 



32j. Accordingly, one can extract A2 from Apy to obtain the strange form 



factors G% + f3G%j from A2, with f3 = rG]^[ / eG'^^^ , if radiative corrections can be neglected. 

To take charge symmetry breaking effect into account, one may simply replace Eq. ( IT2|l 
with 

[1-4: sin^ 6w)GYm - G^e!m ~ - G%^^, (14) 



^E,M 



where G%^^ — | ^G^^^^ - G^'*^^ 



1 / fiu,p 

3 \^E,M 



G 



d,n 
E,M 



and the extraction formula of Eq. 



( IT3l) remains unchanged except G% be re pla ced by G% 



estimated in the constituent quark model 33|, , light-cone meson-baryon model [35| , and 



•,M 



G%M + G%^^. Gg5^f have been 



chiral perturbation theory (xPT) [36] with low-energy constants extracted from resonance 
saturation js^]- 



B. Radiative corrections to the parity-violating ep scattering 

Since the value of A2 in Eq. ( [T3|) is just about a few percent of Ai, it is not possible to 
neglect the electroweak radiative corrections, which is of order 0{aem)i to obtain accurate 
information of the strange form factors of the proton. This is the reason why high precision 



8 



measurements and precise knowledge of the radiative corrections are required to obtain 
reliable extraction of the strange form factors from ep scattering. 

The complete 0{aem) radiative corrections to Apv derive from several different sources 
such as vertex corrections, self-energy insertions of the fermions and gauge bosons, 
7Z mixing, wave function renormalization, two-boson exchange, besides the inelastic 
bremsstrahlung. They have been extensively studied 18N21|. The radiative corrections 



;o Apy have been conventionally taken into account by expressing Apy in following form 



3l| 



Apy{p, k) =v4i + A2 + A3 

Ai = —ap 



1 - 4fi;sin^ 6w) - 



A2 
A, 



ap 



ad -4sin^e 



When the parameters p and k equal one, Eq. (fT5|) reduces to Eq. (j4]), and one recovers the 
tree approximation. The linear combination of the strange form factors, G% + (3G%j, has 
been extracted from A2 in Eq. (ITSl) . In this paper, we will restrict ourself to corrections 
arising from TBE. 



III. THE AMPLITUDES OF TWO-BOSON EXCHANGE DIAGRAMS 



In this section we evaluate the two-boson exchange diagrams in a simple hadronic model 
where the form factors are inserted as regulators and only the nucleon and A(1232) resonance 
intermediate states are included. We present the details of the calculation, including the 
explicit forms of the form factors and the values of parameters employed. As in 25|, |27|, 
28|, we use package FeynCalc |38|] and LoopTools |39| to do the analytical and numerical 
calculations, respectively. 



A. The amplitudes of 2^E and 7Z exchange box diagrams 



Choosing the Feynman gauge and neglecting the electron mass rrie in the numerators, one 
can write down the amplitudes of box diagrams Fig. 1(c) and Fig. 1(d) with the nucleon 
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intermediate states as 

X [(-1 +4sin2eiy] + 75)m(pi^ 



(P4 - ky ~\^ + ie{k- P2Y -Ml + i£ 



X u{p,)Vl{p,, k) ^^''_ ''_^^ Vt{k,p2)u{p2). (16) 

The amplitudes for the cross-box diagrams can be written down similarly. Because the 
amplitudes in Eq. fll6p are infrared divergent, an infinitesimal photon mass A has been 
introduced in the photon propagators to regulate the IR divergence. As explained in |25| . 
in the soft photon limit, the box diagrams and their corresponding bremsstrahlung cross 
section give no correction to Apy- To go beyond the soft photon approximation to estimate 
the corrections to Apy, we calculate the full amplitudes of M*^*^'^) and M^'^'^^ and subtract 
^^soft^ and m]^j'^'* from their respective full amplitude. The interferences between the 
remaining box diagrams and the tree diagrams are then IR safe. 

Similarly, amplitudes for the diagrams Fig. 1(c) and Fig. 1(d) with the A(1232) inter- 
mediate states can be written as follows 



i^TiY (pi + p2 — ky — m1 + is (p4 — kY + is 

-i _ + M^)Pl'hk) . 



(27r)4 ^^^^ ^ + p2 - A;)2 - m2 + M cos Ow ' 

X {{-l+Asm^ew)+l^)u{pi)- 



{P4 - kf + ie{k- P2Y -Ml + ie 
-i[f + M^)P%\k)^ 



X u{j)i)V^^j^^{k,pi - k) — ^2 _ ^2 ^"^g — ^N^Aik, k - P2)u{p2), (17) 
where 

^afi - 9al5 ^ , (J-Oj 

is the spin-3/2 projector. The amplitudes in Eq. ( IT7|) are IR finite because when the four- 
momentum of the photon approaches zero the 7A^A vertices also approach zero. Therefore 
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we do not need to put A in Eq. (fT7|) . The vertex functions T's for A ^ N are defined by 

u{p + q)TT-^NiP,Q>a{p) = -te{N{p + q)\Ji:JA{p)), 

u{p + q)T^^:fj,{p.q)ut{p) = -ig{N{p + q)\J^\Aip)) , (19) 

and similarly vertex functions for — t- A are defined by 

uf{p)K':JAP.QHP'Q) = -^e{A{p)\J:JN{p-q)), 
ut{p)^^Ni^{P.q)u{p-q) = ~ig{A{p)\r,\N{p-q)). (20) 

Note that q's in T'^\l'lf {p, q) and T^l],^^ {p, q) always correspond to the incoming momentum 



of the photon {Z boson), a convention used in 27 1. 



The relations between these vertex functions are 

rA->7v(p. = -7o[r^^A(p. -?)]^7o, rf _^jv(p, = -7o[r|^A(p, -?)]"^7o- (21) 

On the other hand, the following relations 

rA-^7v(p> = 7o[r^_,A(p> 9)]Wo, r|^^(p, g) = 7o[r|_,A(p, ?)]^7o, (22) 

are used in 



271, 



291]. We consider Eq. ( l2Ti) to be the correct one because it can be derived 
from the fact that both of the electromagnetic and neutral weak currents are Hermitian. The 
difference between Eq. (121!) and Eq. (!22|) incurs discrepancies between the results obtained 



m 



28| and [29| as will be discussed later. 



B. Matrix elements of the electromagnetic and neutral weak currents 

between nucleon and A 



Here we discuss the explicit forms of F^^^ and 



magnetic current between N and A is written as 40] 



A^AT- The matrix elements of electro- 



{N{p')\j;r\^{p)) = J^ii(p')[^7iF«(g^)(^?^;^^-p^7"^-7M7>-g + 7M;^g") 
+ 92F^^\q')ip,q''-p-qg^^) 

+ 93F^^\q')/MNiq'ip,r - 9^,^) + q.iq''^ - 7> " g))]75Tjw.(p),(23) 

where q = p' — p and is the third component of the N ^ A isospin transition operator. 
gi are constants and F^[q^ = 0) = 1. One has the following relation between Ge,m,Ci the 
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transition form factors defined by Jones and Scadron [4l|] and gi, g2, gs'- 

3 Mn 



9i 

92 
93 



2Ma + Mn 
3Mn{Ma + SMn 
2 



iGMiO)-GEm. 



Ml- Ml 
M% 



■Ge(0) + ^ 



N 



2Ma + Mn 
Ma + Mn 



-Gc(0) + 



AMI 



M, 



(24) 



2 Ma{Ma + Mn) V Ma-Mn {Ma 
We take G'm(O) = 3.02 [42]. G'ij(O) and GciO) can be inferred from the relations 
Ge{0) = -Gm{0)Rem and Gc{0) = -[AMl/jMl- MI)]Gm{0)Rsai with the experi- 



4^ and R 



SM 



-4.0% 



45|. We thus 



mentally determined values of Rem = —2.5% |43|, 
have Ge{0) = 0.0755 and Gc'(O) = 1.1496 and correspondingly, gi = 1.91, g2 = 2.63, and 
g^ = 1.57. Note that the normalization used in Eq. fl25]) to define the couplings constants 
g'^s differs from that of 



29 



40| where they used Ma instead of M^ everywhere in Eq. 



With this normalization difference taken into account, the corresponding values of g^s used 



m 



29| would be gi = 1.82 and g2 = 2.81, with g^ varied from -0.44 to 1.28. We note that, 



however, since all the current experimental data for Rsm extracted from experiments at low 
Q^'s as small as = 0.060 GeV^ |46| remain negative, we will not consider the possibility 
of a negative value of gs. The difference between the values of gs used in our calculation and 



29| leads to considerable differences in some of the results between these two calculations. 



if the vertex relation of Eq. f l2T]) is employed, as will be discusses in the next section. 
The neutral weak current can be decomposed into isovector and isoscalar parts: 



•^/f — '^vVj^ + I^A^a + isoscalar terms. 



J'^r = + isoscalar terms. 



(25) 



where the superscript "3" refers to the third component in isospin space, ay = (1 — 
2 sin^ 6'vi/)/(2 cos6'vt') and /3a = —1/(2 cos 6'vi/)- The isoscalar part does not contribute to 
— )• A transition. The ZpA'^ vertex contains both the vector and the axial-vector compo- 
nents. The vector part takes the form 

1 



W)\Jlv\^\v)) 



Ml 



+ ~92F'i\q'){v,q''-Vq9l 



+ hFt\q')lMN{q\p,l'' - gli>) + qM'^i' - 7> " q))\l^u^{v). (26) 
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where g'j^s and g[s are related by cji = a/2/3 aygi- Note that the factor yl comes from 
isospin transition operator T3. Thus we have gi = 0.47, ^2 = 0.658 and ^3 = 0.40. 
The axial- vector component of the ZpA^ vertex is given as 28 1 



u{p')[hH^l\q'){g;{p-q)-p,q'^) 



+ h,Ht\q'){gy-p,im<{p). 



+ hH^i\q'){{p-q)ri,-M 



(27) 



where F^^ in Eq. ( 126|) and h''^ in Eq. (l27j) are the vector and axial-vector form factors, 
respectively. In the present investigation, we will assume that, for simplicity, some of them 
separately takes a common form for different couplings, i.e., F^'^ = Fj\{Q'^) and H^^ = 
HaIQ"^)- In addition, both F/^ and if a are normalized to one at = 0. 

Only the coupling constants h[s remain to be determined. They can be obtained from 
the data of uN — > fiA. Many experimental papers on neutrino induced A production adopt 
the notation of Llewellyn-Smith 
axial-current is written as 



(a++(p')I4^'1p(p)) = Mp')f'^^' 



52] where the A^A transition induced by the weak charged 



M, 



N 



Ml 



i{p'-q)9:-qy' 



+ Ct{Q')g 



— 7^p q,4u{p). 



M 



(28) 



N 



The form factors in Eq. fl27|) can be related to the form factors defined in Eq. fl28l) by 
performing a rotation in isospace and assuming the nucleon and A are both on-shell. The 
resulting relations are 



hi 
hs 
hi 



Ma 



V3 



f3cm 



Ml 



Ma 73 



Ml 



+ 



Ml 



(29) 



According to 



47| and 48|], = and hence h^ = 0. If we follow the weak pion production 







50j, then we find (0) 



data of j49| and extrapolate the experimental result to 
—0.8, C^(0) = 2.4, to obtain hi = —0.263, /14 = —0.458. The parameter h2 cannot 
be determined from the weak pion production. According to partial conservation of axial 
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current (PCAC), one has the following relation 



ctiQ') - ^£fQ-.ciiQ% (30) 



where is the pion mass. Hence one obtains ^6^(0) ^ ^C^^iO) ^ 107.7 and the corre- 
sponding value for h2 would be about —360.91. Even with such a large value, we find its 
effect is tiny (< 10^^^) and therefore we simply set /i2 = 0. 

Note that the vertices •jNA and ZNA in Eqs. fl2^ f27} all satisfy the constraints: 

Par^"_^^(p, q') = P/3r^%^(p, q') = 0, (31) 

for any q', to eliminate the coupling of the unphysical spin- 1/2 component of Ratria- 



Schwinger spinor 5l|. The expressions in Eqs. fl26| |27|) have been written in many different 



ways |50l. l52l |53| but only those given here satisfy the above constraints. 

In [29|, different forms of the axial form factors are employed. They obtain the matrix 
elements of J^^ by simply removing 75 from Eq. fl26p and write 



Ma 

It leads to only three form factors instead of four in Eq. fl27|) . The form factor g^iQ"^) of 
Eq. (!32|) is required to have a pole at = 0. However, form factors defined in Eq. ([2] 
are not required to have such poles, and in our opinion, a more appropriate choice. 



C. Nucleon and N ^ A form factors 

So far we have not specified the explicit forms of the nucleon and N ^ A form factors. 
In this article we adopt the following two sets of the form factors. The set A is parametrized 
as follows: 

r'f'P — r^'P In — r^'P i^r — r^'P u, — 1 
^a7^ = (g2 + A2)2' (33) 
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where x = G^/{Q^ = 0), y = G^jfiQ"^ = 0) and z = G^iQ^ = 0). We take Ai = 0.84 GeV 
and A2_= 1.0 GeV from the usual dipole form G^ = 1/(1+QV0.71)2, and Gj, = G^{0)/{1 + 



Q'^Y 121; ISJ], with Q given in unit of GeV , i.e., c = 1, a convention to be used hereafter. 
We determine x, y, z from relations 54 1, G^'^ = p(l — 4/tsin^ (^w)G^^\j — pG% j^.j — pG'^^^ 
and G^ = -(1 + RI=^)Ga + ^R\=^G\ + As at = g point. The quantities G\ and 
As refer to the SU (3) isoscalar octet form factor and the strange quark contribution to the 
nucleon spin, respectively. The p, k and and -R^^" are due to radiative corrections. 

They lead to a; = 0.076 ± 0.00264, y = 2.08 ± 0.00813 - Gl^(O), z = -0.95tUJ^ + As(0). We 
fix X = 0.076 and vary the values of y,z, Ai, and A2 to check the sensitivity of the results 
on the parameters and find little changes. 

The forms of the 7A^A and ZNA are taken to be 



(g2 + A2)2' ) (g2 + AL 



FM') = 77,r^.^ = j^r^-^,. (34) 



Variations of these cutoffs are found not to affect the results significantly as well. 

The form factors set A given in Eqs. ([33D and flMl) do not describe well the existing 
data at large Q"^. For example, the ratio of the proton electric to magnetic form factors 
R = PpGe/Gm has been found to deviate from one at large 2^, while the form factors 
of Eq. (|33|) gives R = 1. Similarly, the N ^ A transition form factors have been measured 
and found to drop faster than Q~'^ at high Q^. More specifically, perturbative QCD predicts 
that at high Q^, the Jones-Scadron form factors scale as follows jsS], 

GuiQ') ~ Q'\ Ge{Q') ~ Q-\ GciQ') ~ g-^ (35) 

such that both Rem and Rsm should approach some constants as — )■ 00. The — > A 
transition form factors given in Eq. (IMjl clearly do not have the correct asymptotic behavior 
at high Q^. We try to take these into account by adding extra factors to both G^'^ and F^^ 
given in Eqs. fl33p and fl34p . This leads to the following more realistic form factors set B, 
withFi^'')(Q2) = FA(Q2), 
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Fitting the data of |2J,|56j gives A3 = 2.0 GeV, A4 = GeV and A5 = 0.5 GeV. Note that 
when one evaluates the effect from the box diagrams with A intermediate states, one still 
needs to specify the choice of the nucleon form factors because one still receives contribution 
from the interference between IjE and TBE box diagrams. Therefore each form factors set 
includes both of nucleon and N ^ A form factors. We will discuss the sensitivity of the 
results with respect to the use of these two different sets of the form factors in Sec IV. 



IV. RESULTS AND DISCUSSIONS 



In this section, we present the results of the corrections of 2'yE and 'yZE to Apy in the 
simple hadronic model described in the previous section. The sensitivity with respect to 
different choices of parameters and form factors will be analyzed in details. The influence of 
the TBE effects on the extracted values of the strange form factors C^^ + f3G%j is discussed 
at the end. 



A. TBE Effects on Apy 



As in 25|, l28(|, we characterize the 2'yE and '~fZE corrections to Apy hj 6 defined as 

Apy{l-f + Z + 27 + 7Z) = Apy(l7 + Z){1 + Sn + Sa), (37) 

where Apy (I7 + Z) denotes the parity- violating asymmetry arising from the interference 
between I7 and Z-boson exchange, i.e.. Figs. 1(a) and 1(b) while Apy{l'y + Z + 27 + 'yZ) 
includes the effects of 2'yE and 'yZE with the nucleon and A(1232) intermediate states. 6n(a) 
represents the contribution from the diagrams with the nucleon (A-resonance) intermediate 
states, respectively. 



1. The TBE corrections from the nucleon intermediate states 

We first present the results of 6]\f as function of e, the contributions of TBE diagrams 
with the nucleon intermediate states, in Fig. |2l The effects of the interferences between 17^' 
and 27ii^ (I7 x 27), and those between IZE and 2'~fE [Z x 27) are represented by the dotted 
and dashed lines, respectively, at four different values, (5^=0.03, 0.1, 1.0, 5.0 GeV^. The 
interferences between 17^' and ^ZE are given by the solid lines. The lines in red correspond 
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FIG. 2: 2jE and jZE corrections to Apy with nucleon intermediate states, as functions of e from 
0.1 to 0.99 at = 0.03, 0.1, 1.0, and 5.0 GeV^, respectively. Dotted and dashed lines denote 
corrections coming only from interferences of I7 x 27, and Z x 27 . The sohd hnes represent 
contributions of 7 x jZ. The hnes in red and black correspond to the results using the form factors 
set A and B, respectively. 

to the results obtained with form factors set A while the black lines are associated with form 
factors set B, as specified in Eqs. ( I33ll36l) in the previous section. We see httle difference 
between red and black curves in Fig. |2l as both form factors sets A and B are of dipole or 
higher order forms. On the contrary, in Fig. [3] one finds that at = 5.0 GeV^ the results 
using the monopole form factors, with cut-offs adjusted accordingly, are much smaller than 



those obtained with sets A and B, as pointed out in [26 1. 

In Fig. 121 we see that both 2'~fE and 'jZE effects strongly depend on and e. The 
magnitude of each contribution has its maximum at e = and decrease to zero when e 
increases. One also sees that I7 x 27 contribution always cancels the Z x 27 contribution 
and hence their sums are always small compared with the size of each contribution, a feature 



also present in the partonic calculation of Ref. [22||. Another interesting fact is that the 
magnitude of SnI^'J x 27) is always larger than 6n{Z x 27). For the I7 x 7Z contribution, 
it decreases as increases and dominates over 6n{'2'jE) at the backward directions when 
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< 0.1 GeV^, but reduces to about the same size as the total contribution at higher 
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FIG. 3: I7 X 27 and Z x 2j corrections to Apy with nucleon intermediate states, as functions of e 
from 0.1 to 0.99 at = 5.0 GeV^. Dotted lines denote the results obtained with monopole type 
form factors. Dashed and solid lines correspond to results obtained with form factors sets A and 
B, respectively. 



2. The TBE corrections from the A(1232) intermediate states 

We continue to present our result of 6^ which arises from the TBE diagrams with the 
A(1232) intermediate states. In Fig. HI we show the 2'-fE and 'jZE corrections to Apy 
by plotting 6a vs. e for both form factors sets A and B. Again, the red and the black lines 
correspond to results obtained with form factors set A and B, respectively. One immediately 
notices that they are very close to each other when < 0.1 GeV^. However, difference 
begins to develop when reaches 1.0 GeV^ at forward angles (e > 0.8). As increases 
further, the difference between red and black lines becomes more pronounced even at small 
e and at = 5.0 GeV^, the discrepancy reaches more that 100% in some cases. The fact 
that 6a is more sensitive than to the details of the form factors indicates that the TBE 
diagrams with the A intermediate states are more strongly dependent on the higher loop 
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FIG. 4: and ^ZE corrections to Apy, with A(1232) intermediate states as functions of e at 
= 0.03, 0.1, 1.0, and 5.0 GeV^. Dotted and dashed lines denote corrections coming only from 
the interference I7 x 27, and Z x 27, respectively. The sohd hnes represent the contribution of 
7 X jZ. The hnes in red and black correspond to results obtained with form factors set A and B, 
respectively. 

momentum than the diagrams with the nucleon intermediate states. 

One further observes that the contributions of I7 x 27 and Z x 27 are neghgible for 

< 0.1 GeV^. As increases, the magnitudes of both contributions increase and become 
comparable in size with 5a (7 x iZ) as reaches 5.0 GeV^. The cancelation between I7 x 27 
and Z X 27 contributions is also seen in Fig. Hlwith the magnitude of I7 x 27 contribution 
larger than that of Z x 27 as in the 6n case. 

The 7 X '-fZ contribution exhibits more complicated and e dependence. At lower 

< O.lGeV^^, it remains small until e reaches between 0.6 ~ 0.8. Then it increases rapidly 
before dropping drastically when e becomes very close to one. For in the region of 
0.1 ~ 1.0 GeV^, 'jZE contribution is fiat and almost zero until e increases past 0.8 and 
becomes small and negative at forward angles. The behavior changes when grows larger 
than 1.0 GeV^, as it decreases monotonically with increasing e, crosses zero at e ~ 0.7, and 
drops rapidly as e reaches 0.9. 
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To sum up, we see that at lower < 0.1 GeV^, 'yZE contribution dominates. When 
reaches 5 GeV^, I7 x 27 effect becomes dominant at backward angles and brings the full 

5a into negative. However, at forward angles the I7 x 27 contribution cancels the sum of 

IZ X 27 and 7 x 7Z and the total 6a becomes negligible. 
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FIG. 5: The short-dashed and sohd Unes correspond to the results using vertex function relations 
Eq. ([22]) and Eq. 1^, respectively. 



3tained with form factors set B. 



Hereafter, we will restrict ourself to results of Sa o 
6a has been calculated independently by [28| and 



29| . The results are different because 



of the following two reasons. The first is that different relations between vertex functions of 
Tn^a and Ta^n are employed. The other arises from employing different values of gs, the 
Coulomb quadrupole excitation strength of — )■ A. 

In general, there are four diagrams associated with the '-/ZE diagram depicted in Fig. 
1(d), two by interchanging the order of the exchanged 7 and Z lines and two others from the 
associated cross-box diagrams. For simplicity, let's just consider only the two box diagrams 
without cross. We denote the amplitude of the diagram with 7 exchanged first by M^z and 
the other with Z exchanged first by Mz-y, both with the A in the intermediate states and 
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only the vector ZNA coupling. We may then write 

3 

s^z = YgjQiCj 



11 

3 

= ^gigjCji- (38) 



Our numerical results for the magnitudes of Cji agree 57| with those obtained in 29|. 



However, with the use of the vertex relation of Eq. ( 122|) . one would obtain 

Cl3 = — C^g, C23 = — C'gg, C3I = —C'^^, C32 = — C'gg, (39) 

where the subscripts 1, 2, and 3 correspond to Ml, E2, and C2 couplings, respectively, and 
Cij = C-j for the rest. On the other hand, there would be no minus signs in Eq. ( I39l) if 
the vertex relation of Eq. ([5T]) is used. Consequently, after summing up 6.yz and 5^^, the 
crossing- couplings of C2 with Ml and E2 terms give no contribution in the calculation of 



29|, while in [28| no cancelation between 6.yz and Sz-y occurs at all. Similar situation also 
takes place with the amplitudes where the vertex ZNA is of axial-vector coupling. The 
resulting discrepancy in the predictions for 6a are shown in Figs. |5] and [6l 

In the upper two figures of Fig. [5l we show the difference arising from using vertex 
function relations of Eqs. ( |2T]) and (1221) for two fixed values of = 0.1 and 3.0 GeV^. 
The sohd and short-dashed lines correspond to the results using Eq. ( 121]) and Eq. (!22|) . At 
low Q^, discrepancy is small. However at higher the difference between two results is 
significant at the forward angles. In the lower two figures of Fig. [5l on the other hand, e 
is fixed at 0.5 and 0.99, respectively, while is varied. We see large difference develop at 
large in both cases. Furthermore, when e is fixed at 0.99, the solid line goes downward 
but the dashed line goes upward, and when reaches 6.0 GeV^ the solid line goes down 
to about -0.04 but the dashed line is almost zero. 



In Fig. in 

as used in 



the dotted and solid lines denote the results obtained with (73 = and (73 = 1.57 



281], respectively. The difference at low cases is very small but becomes 
significant as reaches 3.0 GeV^, especially at the forward angles. It underscores the 



important role played by the Coulomb quadrupole transition in t 



le evaluation of the box 



diagrams with A intermediate states at high and large e. In 29|], (in our convention), 
was varied from —0.44 to 1.28 and the effects of varying are found to be small. It is 
because the use of the vertex relation of Eq. (I22l) leads to cancelations in the cross couplings 
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FIG. 6: The dotted and solid lines denote results obtained with = 0, and g-^ = 1.57, respectively. 

between C2, and Ml and E2 such that the effects of is reduced. Another reason is that 
only cases up to ~ 1.0 GeV^ are explored. 

3. Total effects: sum of 5n and 5a 

Here we compare the behavior of 5n and 5a and present their sum. We see from Fig. [7] 
that both of them are sensitive to and e. At = 0.03 GeV^, 6n- is dominant over 5a 
in the range < e < 0.6 because 6^ is negligible there. As e further increases, 6^ increases 
rapidly at e = 0.6 before dropping at extremely forward angles. These behaviors are in 
sharp contrast with 6n which simply decreases as e increases. The qualitative features of 
the curves at = 0.1 GeV^ remain the same with the one at = 0.03 GeV^. 

When increases to 1.0 GeV^, one sees that 6a is very small and flat while Sn decreases 
monotonously with respect to e. As increases up to 5 GeV^, 6^ becomes negative at 
backward angles but becomes positive as e increases. On the contrary 6n- is always positive. 
We conclude that at small e, 6]\f is dominant but 6^ becomes important as e grows. 

Another way to compare 6n with 5a is to see the evolution of the 6's w.r.t. at fixed 
e as depicted in Fig. [8] for e = 0.5 and e = 0.8. The notation is the same as in Fig. [71 We 
clearly see that for at e = 0.5, 5a is small and of opposite sign to 5Ar, while at larger value 
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of e = 0.8, is always comparable with and becomes dominant at large value of Q^. 

Since 6]y is substantially larger that 6^ in the region e < 0.8, the total effect 6 = 6^ + 6^ 
is very close to 6n- They differ only after e grows larger than ~ 0.8. 



More quantitatively, at = 0.1 GeV^, the full correction, i.e., the combined effect of 6 



N 



and 6a, reaches about 1.75% at backward angle 135° (SAMPLE), about 1.68% at forward 
angle 35° (A4) and about -0.4% at very forward angle 6° (HAPPEX). On the other hand, 
when grows to 1.0 GeV^, the full correction starts from around 1.4% at backward angles 
and decreases to become less than —0.4% at extreme forward angles. 
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FIG. 7: Two-boson exchange corrections with nucleon (dotted) and A (dashed) intermediate states 
to Apv, as functions of e from 0.1 to 0.9 at = 0.03, 0.1, 1.0, and 5.0 GeV^, respectively. The 
solid lines denote their sums 6 = 6n + 



4- Comparison with Marciano-Sirlin approximation 

Here we elucidate the relation between our results with those obtained within MS ap- 
proximation 2o!]. Upon close inspection, the method of MS actually contains three approx- 
imations. The first one is to assume the momentum transfer Q = Pi — = — P2 = ^■ 
Furthermore in the MS approximation the electron mass is neglected and Eiab is taken to 
be zero. This is the second approximation used by MS. Lastly, they take away the Coulomb 
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FIG. 8: Two-boson exchange corrections with an intermediate nucleon (dotted) and A (dashed) 
states to Apv, as functions of from 0.1 to 6 GeV^ at e = 0.5 and 0.8, respectively. The sohd 
hnes denote their sums 6 = 6]\f + Sa- 



interactions because it was argued that its effect has been included in the wave function of 
the bounded electron since they were concerned with the atomic systems. 

Moreover, the MS approximation includes no contribution from the resonance intermedi- 
ate states. Hence we shall compare our results for (5Ar(l7 x jZ) = 6n with values obtained 
in the MS approximation. Since we have already seen that the results do depend somewhat 
on the form factors, we will employ the same form factors used in 20| in order to make the 
comparison more exact. 

We first define some quantities to facilitate the comparison. In the MS approximation, 
the parity asymmetry due to the 7Z diagrams is given by. 



A 



MS,-yZ 
PV 



lab 



(40) 



|MW|2 

where M*^"-' is the \'~)E of Eq. ([3]) and M^^'^^^ the parity- violating part of the '~^ZE amplitude 
evaluated within the MS approximation scheme. The and Eiab dependence of Eq. fHOj) 
arises entirely from M*^") because M^'*^"^ is taken at = and Eiab = 0. 
We further introduce the following quantity. 



OmSW ) ^lab) 



PV,MS-i 
■yZ I 



AO^^iQ^ , ELab) 2Re [M('^)t Mi^)]/ \M(-)\ 



i?e[MWtMW] 



(41) 



On the other hand, the 6n we obtain is given as 



PV,HMi 



6n{Q ,Eiab) 



AOBE [Q2^ELab) 2Re [M Wt M W ] / 1 M(") 1 2 



i?e[MWtMW] 



(42) 
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where M*^''^ is the IZE amphtude of Eq. ([3]) and M^^'^^ the parity-violating part of the 
^ZE amplitude evaluated in our hadronic model with only the nucleon intermediate states 
included, with both dependent on Q"^ and Eiah- The relation between 5^ and 5ms is most 
transparent when = 0, Eiab = because M^^'^^^ is evaluated at this point. In this limit, 
6ms is given [201 as. 



^MS = hlsiQ 

where p^z and K-yz are 
P-yZ = 



0,E, 



lab 



0) = p^z 



4k,^z sin^ 6w 
1-4 sin^ 9 



5ap 



w 



271 



(43) 



2a, 



4 sin^ 6 



n 



w. 



K. 



9 



- 4 sin^ 6 



w 



4 sin^ e 



w 



(44) 



27rsin2 6'vF 

Here K is the asymptotic contribution obtained by carrying out the short-distance expansion 
in a free-field theory. Its value is 8.58 if the onset scale is set to be 1 GeV. On the other hand, 
I^B corresponds to the the long-distance contribution of the 7Z box diagram estimated in 
the Born approximation. Its value has been estimated to be 2.04. 

Hence one obtains Ams^^ow — k) = 1.18% and /S.Ms{high — k)= 4.98%. It was argued in 



291] that the hadronic calculation as done here should correspond to the so-called soft part 



because the form factors used in the hadronic calculation function serves as a regulator and 
the contribution from the higher loop momentum are suppressed. Accordingly, our result for 
6n{Q'^, Eiab) should reproduce Ams{^ow — k) = 1.18% in the proper MS limit as we discuss 
next. 

In Fig. [9] we present our results for 6j\f{Q^, Eiab) by setting = with varying Eiab- The 
full results and the Coulomb contribution are denoted by the solid and short-dashed lines, 
respectively. The difference between the solid and short-dashed curves, represented by the 
long-dashed line, would correspond to the low-k contribution, to be compared with results 



of 



20|. One sees that the long-dashed hue, when Eiab goes to zero, does approach 1.18%, a 



value given in the MS approximation if only low-k contribution is kept. 

In other words, our calculation restores the value given by MS approximation if we follow 
their scheme. On the other hand, it is easy to see that the Coulomb interaction contribution 
is larger as compared with the non-Coulomb part. Furthermore the non-Coulomb contri- 
bution decreases more rapidly as Eiab increases. Note that the calculation in this section 
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FIG. 9: Comparison between our result for 5n = Sn{1^ x 7Z) of Eq. (gS]) at = with 
^MS of Eq. ()43p as obtained within MS approximation. The dashed Hne denotes the contribution 
corresponding to the low-k part of the MS approximation with finite Eiat- The Coulomb part is 
represented by the dotted line and the total contribution is denoted as the solid line. 



is carried out at = 0. We see that the ■yZE contributions is sensitive to Eiab and it is 
necessary to go beyond the MS approximation. 



B. Extraction of the strange form factors 

Here we first examine the effects of the 2'~fE and 'yZE on the values of strange form 

12| . A4 13|, and GO experiments where data have been 



factors extracted from HAPPEX 



taken at forward angles. In SAMPLE experiment, measurements of both elastic ep and 
electron-deuteron (eD) scatterings are combined to extract G^^. However, due to the fact 
that there is no reliable way to estimate the TPE and 'yZE contributions to elastic eD 
scattering, we do not know how to reanalyze SAMPLE data. Naively, one may attempt to 
apply the simple hadronic model here to the deuteron case. But as the deuteron is a loosely 
bound system, treating deuteron in a similar manner as proton is questionable. 
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1. The formulation of the extraction of the strange form factors 



All the existing analyses ISMlTf extract strange from factors from Eq. ( IT5|) where the 



electroweak radiative corrections are included in the parameters of p and k in the expression. 



The latest PDG values [58[ for p and k are p = 0.9876 and k, = 1.0026. They deviate from one 
because higher-order contributions like vertex corrections, corrections to the propagators, 
and TBE effects are taken into account. Since we have explicitly calculated the effect of 
TBE effects in this study, we should replace the contribution of TBE to the above-mentioned 
p and K as estimated by MS, with our results to avoid double counting. Namely, one should 
then subtract Ap = p^z and An = k,^z from p and n and use p' = p — Ap and k' = n — An 
in Eq. ( fTSi) instead. 

As explained earlier, p^z and k^z in Eq. (jUj) actually consist of two contributions, 
namely, the high-k and low-k parts and our results correspond to the low-k part only. We 
should then only take away the soft loop momentum contribution, which is associated with 
the term, and define Ap, An as follows: 

Ap = -^(1 -4sin2^vK) ■ -^B = -0.73 x 10"^ 

TT 5 

Q, /g \ 4 

An = ^f— - - 4 sin^ (1 - 4 sin^ Ow) ■ -^b = -1.03 x 10^1 (45) 

27ism Ow \4 J 5 

Consequently, we set the experimental parity asymmetry A^py^'' as follows: 

^(f^) ^ Apy(l7 + Z + 27 + 7^), 

= Apv{p\k'){1 + 6). (46) 

With the value we obtain for 6, we can determine Apy{p', k') and extract the strange form 
factors from the resultant A2 of Eq. ( fT3l) . 
Furthermore, we introduce 



G% + PGl, = iG% + PGl,)il + 6g), (47) 

to quantify the effects of the 2'yE and 'yZE to the extracted values of G|; + l3G\j, where 
G% + [iG%j and G% + [iG^j are extracted from Apv{p, n) and Apv{p', n'), respectively. From 
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we get, 



- A, 



A Exp 



(Aap sin^ 9w 
- A, 



( - A. 
\A%^-A^ 




Ap 
P 



(49) 



where Aq = Ai{p, k) + A^. Note that the values of r/i, r]2, and 773 all depend on the values of 
the inputs of the nucleon form factors such as G^'^^, G'^^^i ^a- -^^ ^ result, the value 
of 6g also depends on those inputs. 

201 within 



We further define 60, the corresponding value of 6 as would be obtained in 
Q = approximation scheme such that Sq = if S = 6q. In other words, difference between 
S as we obtain and 60, represents the possible Q^-dependence neglected in the estimation of 
13], such that 6g vanishes when 6 = 60. Explicitly the value of 60 is given by 



5o = -^A.-!^^. 



(50) 



Obviously the value of 60 depends on the inputs of the proton and neutron electromagnetic 
form factors as well. 



2. Extraction of the strange form factors at HAPPEX, A4, and GO experiments 



At forward angles, A3 in Eq. (fT5|) is negligible because both e' = a/t(1 — r)(l — e^) ^ 1 
and 1 — 4sin^^vy ^1- It offers some advantages that the strange form factors can then 
be determined more accurately. This is why many experiments, like HAPPEX, A4, and GO 
are carried out at very forward angles. In Table HI w e present our results for 6n, 5a, their 
sum 6, besides 6q and 6g, for HAPPEX 12|, A4 13], and GO jl^ experiments. They are 
obtained with the use of form factors set B. We also list the values of Gg = G% + f3G%j. 



For the GO experiments, only measurements of Apy are given in [Ij] and the corresponding 
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values of Gg listed in Table |T] are extracted by us with the use of the nucleon electromagnetic 
factors parametrized in [59]. The resultant change in the values of Gg after TBE effects are 
properly taken into account, i.e., AGs = {G% + /3G%j) — {G% + /3G%[), are given in the last 
column of Table [H 



Exp 




e 


SNi%) 


<5a(%) 


6{%) 


So{%) 


5g(%) 


G.(10-2) 


AG,(10-2) 


HAPPEX 


0.477 


0.974 


0.18 


-0.27 


-0.09 


0.20 


-2.54 


1.4 


-0.04 


HAPPEX 


0.109 


0.994 


0.21 


-0.80 


-0.58 


0.51 


-20.63 


0.7 


-0.14 


GO 


0.122 


0.9930 


0.21 


-0.72 


-0.51 


0.61 


-3.63 


3.9 


-0.14 


GO 


0.128 


0.9926 


0.21 


-0.70 


-0.49 


1.05 


-1.28 


9.2 


-0.12 


GO 


0.136 


0.9921 


0.21 


-0.67 


-0.44 


0.81 


-1.60 


7.7 


-0.12 


GO 


0.144 


0.9916 


0.20 


-0.64 


-0.41 


0.38 


14.14 


-1.1 


-0.16 


GO 


0.153 


0.9911 


0.20 


-0.61 


-0.39 


0.51 


-3.50 


3.8 


-0.13 


GO 


0.164 


0.9904 


0.20 


-0.58 


-0.36 


0.41 


-9.18 


1.5 


-0.14 


GO 


0.177 


0.9896 


0.20 


-0.55 


-0.32 


0.31 


6.19 


-2.3 


-0.14 


GO 


0.192 


0.9886 


0.19 


-0.52 


-0.29 


0.35 


-16.05 


0.8 


-0.12 


GO 


0.210 


0.9875 


0.19 


-0.48 


-0.29 


0.30 


48.25 


-0.3 


-0.14 


GO 


0.232 


0.9860 


0.19 


-0.44 


-0.25 


0.30 


-20.25 


0.6 


-0.12 


GO 


0.262 


0.9840 


0.19 


-0.40 


-0.21 


0.35 


-2.26 


4.6 


-0.10 


GO 


0.299 


0.9814 


0.19 


-0.36 


-0.17 


0.26 


-8.68 


1.2 


-0.10 


GO 


0.344 


0.9783 


0.19 


-0.32 


-0.13 


0.28 


-1.99 


4.4 


-0.09 


GO 


0.411 


0.9735 


0.19 


-0.27 


-0.08 


0.27 


-1.18 


6.4 


-0.08 


GO 


0.511 


0.9657 


0.20 


-0.23 


-0.03 


0.19 


-2.10 


2.8 


-0.06 


GO 


0.628 


0.9580 


0.21 


-0.20 


0.01 


0.20 


-0.71 


6.8 


-0.05 


GO 


0.786 


0.9413 


0.22 


-0.18 


0.04 


0.15 


-0.81 


3.9 


-0.03 


GO 


0.997 


0.9197 


0.25 


-0.18 


0.07 


0.15 


-0.32 


7.6 


-0.02 


A4 


0.108 


0.83 


1.07 


0.53 


1.60 


0.61 


2.00 


7.1 


0.14 


A4 


0.23 


0.83 


0.66 


0.14 


0.80 


0.29 


2.85 


3.9 


0.11 



121, GO 



]J] , and A4 



TABLE I: The values oi 6n,6a, and their sum 6 for the HAPPEX 
We give the values of 6o, 6g, Gg, and AGg obtained with = 1.57, for those data points. 



13( 1 data. 
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All experimental data included in Table |T] were obtained in the near forward directions 
with e > 0.8. More specifically, the HAPPEX and GO data were taken at extremely forward 
angles with e > 0.92. It is seen from Table [Ithat in this region there is a cancelation between 
6n and 6a as they are of opposite sign. The magnitude of 6a is always larger than 6n at 
lower < 0.7 GeV^. When increases past 0.7 GeV^, 6n overtakes 6a and the sum 6 
becomes positive. On the other hand, in the kinematical regions of A4 data, both 6^ and 
6a are positive such that the sum 6 is also positive. 

Furthermore, one notices that 6q is in general larger than 6^. It implies that MS approx- 
imation overestimates the TBE contribution, besides neglecting the strong dependence. 

The values of 6g presented in Table [T] are considerably smaller than what we reported 



m 



. It can be understood from the following reasons. First, the values of 6's listed are 



already different from before since they are obtained with different nucleon and N ^ A 



form factors. Previously in [25|, |28[, the nucleon form factors used are of monopole type 
while the A^A transition form factors are taken to be of dipole form. In addition, p^z and 
K.yz obtained here are now used to replace only the soft part contribution evaluated by MS, 
i.e., the Ap and An of Eq. ( H5l) . as emphasized in 29|. For example, for the HAPPEX data 



at Q'^ = 0.109 GeV^ and e = 0.994, the value of 6g in [28| is given as -75.23%. However 
if we use the value of 5 = —0.58% in Table [T] instead of the previous 6 = —1.19%, then 6g 



is reduced to —60.52%. I 
hadronic contribution in 



we further use the value of n^z = —1.03 x 10 identified as the 



20l | instead of the value of —5.33 x 10 then value of 6g becomes 



—14.99%. Lastly, using the value of Ap = —0.73 x 10^^ given in Eq. ( H5l) . instead of the 
value of —3.72 x 10~^, produces small change and leads to the final value of 6g = —20.63% 
as given in Table I. 

In general, the values of 6g are smaller than 10% and are mostly negative with the 
exception of backward data of A4. For HAPPEX data at = 0.109 GeV^ and GO data at 
Q'^ = 0.144, 0.192, 0.210 and 0.232 GeV^, the magnitudes of 6g are large and range between 
—20.95% to 63.73%. The magnitudes of 6g seem to behave irregularly. However if one 
computes the AGs = {G% + P&m) — + f3G\j), the values of AGs are relatively stable 
with typical size of —(0.1 ~ 0.2) x 10~^. It is because those with large values of 6g have 
small values of Gg- 

Lastly, to illustrate the sensitivity of the corrections to the extracted strange form factors, 
with respect to the possible experimental uncertainties in the extracted value of Rsm and 
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the resulting Coulomb quardrupole excitation strength of the A(1232), we give in Table II 
our results for 6^, S 60, and 6g, obtained with = for some of the HAPPEX, A4, and GO 
data. Comparison between Tables I and II shows that the variations in the final corrections 
to the extracted values of the strange form factors, when gs changes from to 1.57, amount 
to about 20%. 



Exp 




e 


Sa{%) 


6{%) 


Sg{%) 


AG,(10-2) 


HAPPEX 


0.477 


0.974 


-0.30 


-0.13 


-2.81 


-0.04 


HAPPEX 


0.109 


0.994 


-0.94 


-0.73 


-23.36 


-0.16 


GO 


0.128 


0.9926 


-0.82 


-0.61 


-1.39 


-0.13 


GO 


0.144 


0.9916 


-0.75 


-0.55 


16.05 


-0.18 


GO 


0.164 


0.9904 


-0.68 


-0.48 


-10.32 


-0.15 


GO 


0.210 


0.9875 


-0.56 


-0.37 


54.38 


-0.16 


A4 


0.108 


0.83 


0.58 


1.65 


2.11 


0.15 


A4 


0.23 


0.83 


0.17 


0.83 


3.03 


0.12 



TABLE II: The values of 6a, 6, 6g, and AGs obtained with 93 = for some of the HAPPEX, A4, 
and GO data. 



V. SUMMARY 



In summary, we present the details of our calculation [25|, |28| of the two-boson exchange 
effects in the parity-violating ep scattering within a simple hadronic model with both the 
nucleon and A(1232) resonance intermediate states included. We examine the sensitivity of 
the results with respect to the form factors. We find that the nucleon contribution 6n does 
show mild sensitivity to the form factors depending on whether monopole or dipole form 
factors are used. However, little difference is found between results obtained with a purely 
dipole form factors set A and another more realistic form factors set B which differs from 
set A only at higher Q^. For the A contribution 6 a, however, predictions obtained with the 
use of form factors sets A and B do exhibit substantial difference at high Q^. 

In addition, we compare our calculation [28] for 6a with a recent calculation of Ref. 
291] where different relations relating vertex functions of T^^a and T^^jy are employed. 
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Considerable discrepancy shows up at Q"^ > 3.0 GeV^ and e > 0.5, when the Coulomb 
quardrupole excitation (C2) strength of the A, is nonvanishing. Accordingly, if one takes 

= 1.57, a value determined from the recent pion electroproduction data [4^ is used, our 
results for (5a differ significantly with those given in {29]. 

Furthermore, we clarify the relation between our results and the well-known results of 
the 'jZE effects given by Marciano and Sirlin (MS). We explicitly demonstrate that our 
calculation, with only nucleon intermediate states included, restores the values given by 
MS as long as we follow their scheme to set Q = 0, Eiab = 0, and remove the Coulomb 
interaction. 

We find that both the nucleon contribution 6n and A contribution depend on both 
and e. 6n is always positive and decreases with increasing e. On the contrary, A contribution 
5a exhibits stronger dependence on both and e. In general, 6^ dominates over 5a except 
at extreme forward angles. The sum 5 = + (5a is then positive for e < 0.95 and turn 
negative after then. 

We also present our result of the correction to the extracted values of the strange form 
factors G% + PG%j from the HAPPEX, A4, and GO data at forward angles. Comparing 
with the previous result 25|, |28|, the updated values are reduced. However, the modification 



incurred in going beyond the MS approximation is still significant (up to ~ 60%) for some 
data. In addition, the sensitivity of the correction to the extracted C^^ + PG'l.j values with 
respected to the experimental uncertainty in the determination of Rsm is found to give rise 
to about 20% variations when Rsm changes from to —4.0%, or equivalently (73 = ~ 1.57. 

As we find significant contribution from TBE with A excitation in the extreme forward 
direction, where many of the current experiments are performed, question of the inclusion of 
higher resonances comes up naturally. Naively, one would expect that A(1232) would give 
the largest contribution since it is the most prominent resonance at low energies. Higher 
resonances would be suppressed because of their larger masses. However, only explicit 
calculation can answer this question. Recent dispersion relation calculation of the 'jZE 
correction to Qw |6Q|] could be used to clarify this question in the exact forward scattering. 
However, our results indicate that 6 depends sensitively with at low momentum transfer 
so whether dispersion relation method of [60i] can be extended to investigate the TBE cor- 



rection to strange form 
use of GPD as done in 



actors remains to be further explored. Study of TBE effect with the 



221] and 23|] for TPE effects, will also be very helpful in this regard. 
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